Surface defreezing of glasses 
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A glass surface may still flow below the bulk glass transition temperature, where the underlying 
bulk is frozen. Assuming the existence at T = To of a bulk thermodynamical glass transition, we 
show that the glass- vapor interface is generally wetted by a liquid layer of thickness I ~ — ln(Tb — T) 
when T — > T ~. Contrary to standard surface melting of crystals however, the integrated value of 
the diffusivity across the interface remains finite for T — » T ~ . Difference in shape induced by bulk 
and by surface flow is discussed as a possible means of experimental detection of surface defreezing. 
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Glasses embody the paradigm of broken ergodicity in 
condensed classical systems. While most approaches to 
glasses address an infinite, homogeneous system, there 
is a clear scope for extension to inhomogeneous situa- 
tions, as they may be variously realized in real life. The 
most common -and conceptually the simplest- kind of 
inhomogeneity is represented by the surface. The prob- 
lem we wish to address here is the state of a free glass 
surface, in the neighborhood of the bulk glass transition 
temperature To. 

In crystals near the bulk melting temperature Tm, sur- 
face melting is well documented both experimentally and 
theoretically. As T — > Tm along the solid-vapor bulk 
coexistence line, most crystalline faces of a majority of 
substances develop a microscopic liquid film which spon- 
taneously wets, in full thermal equilibrium, the solid- 
vapor interface The thickness of the film diverges 
as T — > Tm, explaining among other things why crys- 
tals with free surfaces cannot sustain overheating. The 
surface acts as a 'defect' where the liquid phase micro- 
scopically nucleates. 

Due in part to the lack of a comparable understand- 
ing of the properties of even bulk glasses, the possibility 
that similar phenomena might take place at the surface 
of glassy materials has been paid very little attention so 
far, in spite of its potential importance, both conceptual 
and practical. Conceptually, the purely dynamical arrest 
typical of glasses provides a neater example of broken 
ergodicity than that of the solid, where the crystalline 
order parameter represents an additional complication. 
Practically, the possible flow of surfaces can be expected 
to play a role in measurable properties of glasses, such 
as friction, surface flow under intense acceleration, inter- 
facial contact, and other phenomena. In this letter, we 
base on a thermodynamic glass transition theory, of the 
type recently considered within the glass community |jj , 
our first attack to the glass surface problem. 

Given a thermodynamic formulation with a well de- 
fined free energy versus some order parameter, such as 
the atomic density n(r), or the energy density e(r), the 
natural approach to try first is a Landau theory [Q . For 
crystals, it provides what certainly is the simplest micro- 



scopic theory of surface melting, summarized as follows. 
Call f(n) the bulk free energy density, a function of the 
(uniform) atomic density n. The global free energy cost 
of a solid- vapor interface from the solid at x = —00, to 
the vapor at x — +00 is 



F[n(x)} 



[f(n(x)) + (J/2)(dn(x)/dx) 2 ] dx (1) 



Strictly short range forces are assumed, and the gradient 
term (J > 0) accounts as usual for reluctance of the or- 
der parameter n against spatial change, assumed to be 
sufficiently slow. Along the solid-vapor coexistence line, 
and just below the triple point Tm, f(n) will exhibit, be- 
sides the two identically deep (solid and vapor) minima 
f(n s ) — f(n v ) at n = n s and n v respectively, a third, 
shallower minimum f(n) — f(n s ) + A at the intermedi- 
ate liquid density n l , in which A = £q(Jm — T), where 
£0 > is a constant. 

Minimization of the free energy (Q) is formally identi- 
cal, via the dictionary n — > z, x — > t, to minimization of 
the action of a classical point particle of coordinate z and 
mass J, moving in a potential V{z) = —f(z), with energy 
exactly equal to —f(n s ) = —f(n v ). The resulting equi- 
librium interface profile n(x) exhibits, between the solid 
and the vapor, an intermediate liquid film. Calling 7 the 
curvature of f(n) at the liquid minimum, the thickness 
of the liquid film is given by I ~ — \J J/7 ln(A), which 
diverges as T approaches Tm from below, demonstrating 
surface melting (in real crystals, the logarithmic diver- 
gence usually turns to a power law, due to long range 
forces, not included in (|l|) [jl|,||). 

To apply arguments similar to these to a glass-vapor 
interface we must have, as a starting point, a thermo- 
dynamical description for the properties of bulk glasses. 
For that, a first necessary assumption is that crystalline 
states of the system, even if lower in free energy than any 
of the glass configurations, can be ignored. Experimen- 
tally, crystallization of glasses can be kinetically avoided 
by a sufficiently rapid cooling in absence of crystalline 
germs. With this assumption we can consistently speak 
of the glass transition as an equilibrium phenomenon, 
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even though it actually occurs on a metastable branch of 
the phase diagram. 

We will thus suppose that there is a thermodynamic 
glass transition for the bulk system. The configurational 
entropy of the system s c as a function of the enthalpy 
h (to be used instead of the internal energy e since we 
work at constant pressure) should vanish for h lower than 
some ft-o, and increase linearly for h greater than hq. The 
glass transition is in this scheme a second order, mean 
field transition [Q. The critical temperature To is given 
by T^ 1 = ds c /dh\ hg . For T < T the system freezes 
in the configurational ground state, and thermodynami- 
cal variables such as specific heat give information only 
about the vibrational structure of the valley around the 
ground state. For T > T a configurational contribu- 
tion arises, that counts the number of different valleys 
that the system is able to sample. The specific heat has 
a finite jump at T = To, sometimes referred to as the 
Kauzman temperature |^,^] . 

A further connection of this thermodynamic picture 
with dynamical properties at equilibrium is given by 
the phenomenological Adam-Gibbs formula 0, that re- 
lates the configurational entropy to dynamical variables 
such as the viscosity, or the diffusivity D, in the form 

D = £>oexp ^— yf - ) wnere Dq > and A > are con- 
stants. This formula can be made heuristically plausible 
but is not rigorous, and should be considered just 
a useful working hypothesis. As T — > To, s c and thus 
D vanish, reflecting the impossibility for the system to 
jump between valleys, for only one is thermodynamically 
favored. When applied to the previous description of a 
bulk glassy phase, the Adam-Gibbs formula is consistent 
with a Vogel-Fulcher-Tammann type temperature depen- 
dence of the diffusivity, which is experimentally well ver- 
ified for many different glass formers [9 . 

Next, we will assume our system to possess in addition 
to the glass-liquid transition, a liquid-vapor transition. 
At each temperature below the liquid-vapor critical tem- 
perature the free energy f{h) will have two minima, one 
corresponding to the vapor and the other to the con- 
densed phase, the same minimum comprising both liquid 
and glass because, as Fig. [j] indicates, the glass-liquid 
transition is second order. At the very minimum the sys- 
tem will be liquid or glassy depending on whether T is 
lower or larger than To. The second order nature of the 
glass-liquid transition reflects in a jump in the second 
derivative of f(h) at some critical value of the enthalpy 
h c (h c may be different from ho due to vibrational con- 
tributions). On the phase boundary between the vapor 
and the dense phase the two minima of the free energy 
are degenerate, f{h\) = f{h 2 ). 

We can now model the glass-vapor interface at bulk 
coexistence, where the bulk enthalpies of the two coex- 
isting phases are hi and hi. We will suppose, as in the 
crystal- vapor interface, that the system can be assigned a 



well defined, slowly varying enthalpy density h(x) at each 
position x across the interface. Moreover, for x — > — oo 
(+oo), h — ► h\ {h?), forcing the existence of the interface 
through the boundary conditions. In the end, we will de- 
rive the interface profile h{x) by minimizing a free energy 
functional qualitatively similar to that in (Q), but using 
enthalpy as the order parameter instead of density, and 
with only two minima instead of three. 
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FIG. 1. (a) Generic pressure-temperature phase diagram 
for a system with glass, liquid, and vapor phases (a possible 
crystalline phase is not indicated) . (b) Free energy density as 
a function of the enthalpy at two coexistence points on the 
first order line of (a). The coexistence is between liquid and 
vapor in (1), and between glass and vapor in (2). 




FIG. 2. Evolution of the interface enthalpy density from 
the glass at the left to the vapor at the right (T < To). The 
glass-vapor interface is coated by a liquid film of thickness 
I ~ Xo — x. 

A schematic representation of the interface enthalpy 
profile h{x) in the case T < To is given in Fig. |^. This 
function interpolates between the two minima hi and hi 
of Fig. |l|(b). We expect a weak singularity of h{x) where 
h = h c , the point in which f{h) has a jump in the second 
derivative. The corresponding point x marks the border 
between the glass at the left, an the liquid at the right. 
The liquid transforms into the vapor phase roughly at 
some x = xo, where the maximum of f{h) is overcome, 
and the rate of change of h{x) should be maximum. The 
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difference I = xq — x can be identified with the thickness 
of the liquid layer that wets the glass-vapor interface. 
We are interested in the behavior of the system around 
T = Tq. Exactly at this point h c coincides with the 
minimum of /, and for slightly different values it may be 
supposed to have a linear dependence on Tq — T, namely 
h e -h 1 = e(T -T). 

To calculate h{x) and from that I, we proceed in the 
following way. As already mentioned, the free energy 
f(h) is smooth except in the point h c , where it has a 
jump in the second derivative. We will approximate f(h) 
with its quadratic expansion close to h c 



f{h)=ih 2 /2 

f{h)= 1 {h-h c {l-ih)) 2 l2 + .f a 



h < h c 

h> h c (2) 



where we have taken h± — for simplicity. 7 and 7' ( 
7 > 7') are constants related to the bulk properties of 
the material, and /q is a constant ensuring continuity of 
/ at h c . The Euler equation satisfied by the equilibrium 
profile h{x) is, from (fil) 



df(h)/dh = Jd 2 h/dx 2 



(3) 



For f(h) given by (|^), the solution can be written in term 
of exponentials, and a direct calculation shows that 



h(x) = Ae^~ Jx 

h(x) =u+ (h m - lOeV^ 7 * 



h < h c 

h > h c (4) 



where 

u = h c (l ~ Hi) and A = h c [{h m - u)/(h c - u)}^^ 
are values that depend on temperature through h c . Here 
we have set xq = 0, assuming that the barrier of f(h) 
is parabolic up to the maximum, that occurs precisely 
at h m . Then the solution (^) describes the interface in 
the condensed region. Deviations are expected to oc- 
cur around the liquid-vapor transition. The position x 
of the liquid-glass interface, namely the point at which 
h(x) = h c can be written as 

x = -y/Jfrln [(h m - h c (l - ih))/{hci/l)] (5) 

This value diverges when T — > Tq ', indicating the pres- 
ence of a diverging liquid layer of thickness 



e(r -r)V 



(6) 



The spatial scale for the thickness of the liquid layer is 
thus given by yj J/7 , exactly as for a liquid layer at 
the crystal-vapor interface. This underlines a common 
physical origin, namely a lower surface free energy of the 
liquid relative to the solid. 

The logarithmic increase of the liquid layer thickness 
close to To is at first sight also the same as that of the 
liquid layer in short range crystal-vapor surface melting. 



However the conceptual origin of the logarithmic behav- 
ior is different here, as it relates to a second order liquid- 
solid transition, as opposed to a first order one in the 
crystal case. In practice too, glass surface melting should 
be quite different. In crystals, the density jump associ- 
ated with melting gives inevitably rise to different optical 
conductivities of solid and liquid, leading to a nonzero 
Hamaker constant H and to long range dispersion forces 
~ H/l 2 A positive Hamaker constant will thus en- 
hance surface melting, and generally transform the loga- 
rithmic into a power law growth, / ~ (Tm — T)~ 1//3 
A negative Hamaker constant will instead suppress, or 
block, surface melting ||. So will, for a different rea- 
son, commensuration of surface layering with the spacing 
of crystal planes |10|. In the glass there are no crystal 
planes to block surface melting. Moreover, the bulk den- 
sity and optical conductivity of glass and liquid are not 
expected to differ discontinuously and the Hamaker con- 
stant should basically vanish. In conclusion, in glass sur- 
face melting the logarithmic film growth behavior should 
be more robust. 

We can also calculate the diffusivity profile D(x). The 
configurational entropy vanishes linearly with h close to 
the liquid-glass transition, namely 



s c (x) = a{h(x) - h c ) = 



a 




h m -h c (l - 


- he- 








ir. 



(7) 



where a is some constant. Using the Adam-Gibbs for- 
mula, we obtain 



ln(D(x)/D ) 



A 
Ta 



e(T - T) 



, (8) 



valid for T close to To, and x close to the liquid-glass 
transition value i. For larger values of x — x, the non- 
quadratic nature of the free energy functional (partic- 
ularly the maximum between the vapor and condensed 
phases) should be taken into account properly. If T < Tq, 
the previous expression is valid for x > x, and D is for 
x < x. If T > To, i.e., if the bulk system is in the liquid 
phase, then (g) is valid for all x. In this case D{— 00) 
is different from zero, and (^) predicts an enhancement 
of the diffusivity at the surface compared to that of the 
bulk. We have plotted the behavior of D(x) in Fig. || for 
different temperatures. It should be noted that the inte- 
grated diffusivity D = J D(x)dx tends to a finite value 

D - D o yj/7exp ( To ~fa m ) when T — > T ~, This is at 
variance with what happens for melting of crystal sur- 
faces, where D ~ I. Hence surface melting is in a sense 
much weaker on the glass surface than on the crystal 
surface. 
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FIG. 3. Diffusivity across the glass-vapor interface for dif- 
ferent values of A/h m = eo(Tm — T)/h m , indicated on the 
curves (the particular value y' /f = 0.5 was used). The dense 
phase (glass for A > 0, liquid for A < 0) is at the left, while 
for a;>0 the system is in the vapor phase. The curve for the 
bulk glass transition temperature T = To is shown as a dashed 
line. Note the finite diffusivity in a "melted" surface film for 
T < T , and also the residually higher surface diffusivity just 
above T — To. 

Experimental detection of the liquid layer predicted by 
theory to exist at the glass- vapor interface should be pos- 
sible. One method might be to look for "slumping" (thin- 
ning of one end and thickening of the other end) of ini- 
tially shape-controlled glass samples in a centrifuge. The 
shapes expected for bulk-flow-induced and for surface- 
flow-induced slumping are different. Generally speaking 
bulk flow under acceleration should preserve sharp edges 
and flat profiles, transforming e.g., a rectangular shape 
to a trapezium H]. Surface flow under parallel acceler- 
ation will modify a given profile y(x) in time according 
to an equation of the type 113] 



dy _ n 



y 



(l+y*)3/2 



(9) 



where y' = dy/dx, and C is roughly proportional to D 
and to acceleration g. This sort of nonlinear heat conduc- 
tion evolution will not generally preserve sharp edges and 
flat profiles. Given for example an initial step function 
y(x, t = 0) — yo6(x) , it will evolve in the following man- 
ner. First, both edges, the upper one yo < y < yi(t), and 
the lower one y2(t) < y < will become smeared, while a 
central window y 2 (^) < y < Vi{t) °f the face will remain 
flat. As the spatial extension of the smeared corners in- 
creases with time, the window will gradually shrink, and 
after a critical time t c it will close, eliminating all traces 
of a flat face in the slumped profile. The time to required 
for the surface smearing front to advance by yo is roughly 
given by to ~ yoVo/(loPg), where rjo is the viscosity of the 



superficial defrozen film, i.e., a value typical of the liq- 
uid, p is the density, and Iq ~ yj J/7 is the film thickness. 
The possibility of observing this effect (macroscopically, 
by reflectivity, and even better microscopically, by some 
surface topographic technique) seems quite plausible. 

In surface sensitive calorimetry, alternatively, a pro- 
gressive defreezing of the surface film could yield a de- 
tectable mark. Surface melting of crystals yields a 
singularity in the specific heat of the form ~ (Tm — 
T)-i-i/(^-3) when T approaches T M from below, (if 
long-range forces decay as l/r v ). In the glass, inte- 
grating the enthalpy profile (|4|) across the interface we 
obtain a singularity in the specific heat of the form 
- T) when T -> IT . This is again 
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a divergence, as in crystal surface melting, but now a 
much weaker one. In compound or polar glasses, the 
surface static dielectric permittivity could also be a use- 
ful tool. Microscopically, finally, liquid-like surface flow 
could be detectable by surface-specific techniques, for ex- 
ample thermal atom scattering. 

Glass surfaces can also be very effectively simulated. 
A simulation study of PdsoSi2o published some time ago 
does provide a first qualitative suggestion for surface de- 
freezing in a realistic glass forming system Jl3|] . The glass 
surface roughness drop demonstrated in that work upon 
cooling could now be further elaborated upon in terms 
of suppression of capillary fluctuations due to confine- 
ment of the liquid film. Our present mean-field theory 
does not include capillary fluctuations. Their effect on 
the liquid film thickness is known to be at most marginal 
in three dimensions in crystal surface melting |l4|], but 
will require a specific treatment in the glass case. Exper- 
imental evidence of surprisingly flat glass surfaces pl| do 
exist and might be related, although the effect of grav- 
ity is not discounted. This is a line of research that will 
deserve further theoretical and experimental effort. 
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